Magnetic Phase Transitions in One-dimensional Strongly Attractive 
Three-Component Ultracold Fermions 
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We investigate the nature of trions, pairing and quantum phase transitions in one-dimensional 
strongly attractive three-component ultracold fermions in external fields. Exact results for the 
groundstate energy, critical fields, magnetization and phase diagrams are obtained analytically from 
the Bethe ansatz solutions. Driven by Zeeman splitting, the system shows exotic phases of trions, 
bound pairs, a normal Fermi liquid and four mixtures of these states. Particularly, a smooth phase 
transition from a trionic phase into a pairing phase occurs as the highest hyperfine level separates 
from the two lower energy levels. In contrast, there is a smooth phase transition from the trionic 
phase into a normal Fermi liquid as the lowest level separates from the two higher levels. 
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There is considerable interest in three-component ul- 
tracold fermions [l], 0, 0, 0] . Atomic Fermi gases with 
internal degrees of freedom are tunable interacting many- 
body systems featuring novel and subtle quantum phase 
transitions. Two-component Fermi gases of ultracold 
atoms with population imbalance have been experimen- 
tally observed to undergo a quantum phase transition 
between the normal and superfluid states Q . The bound 
pairs form a Bardeen-Cooper-Schrieffer (BCS) super- 
fluid, while the unpaired fermions remain as a separated 
normal Fermi-Liquid (FL). These exotic phases have re- 
vived interest in the one-dimensional (ID) integrable 
model of two-component fermions, which captures the 
physics involved in quantum phase transitions and mag- 
netic ordering [f| 0, 0, 0] . 

Three-component fermions reveal more exotic features 
0, 0, 0, II EI El E3 El- The scattering lengths be- 
tween atoms in different low sublevels are again tunable 
via Feshbach resonances [3, EH EH • As a consequence, 
BCS pairing can be favored by anisotropies in three dif- 
ferent ways: specifically, atoms in three low sublevels 
denoted by |1), 2} and |3) can form the three possible 
pairs |1) + |2), \2) + |3) and |1) + |3) 14]. One may also 
have three internal spin states exhibiting SU(3) symme- 
try via tuning three scattering lengths close to each other 
17]. Significantly, strongly attractive three-component 
atomic fermions can form spin-neutral three-body bound 
states called trions. Thus, a phase transition is ex- 
pected to occur between pairing superfluid and trionic 
states P, 0,0, Ell- 
in this Letter, we consider ID three-component ul- 
tracold fermions with ^-function interaction in external 
magnetic fields. Although this model was solved long ago 
by the Bethe ansatz (B A) [l8|, El| , its physics is far from 



being thoroughly understood. Here we study the precise 
nature of trions and pairing in this model, and calculate 
critical fields and full phase diagrams by solving the BA 
equations and related dressed energy equations. Our an- 
alytical results for magnetism and magnetic-field-driven 
quantum phase transitions in attractive fermions should 
provide benchmarks for experiments with ultracold Fermi 
atoms with multiple internal states. 

The model. The Hamiltonian 18| we consider 
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describes N fermions of mass m and spin-independent s- 
wave scattering lengths, which can occupy three possible 
hyperfine levels (|1), \2) and |3)) and are constrained to a 
line of length L with periodic boundary conditions. The 
last term denotes the Zeeman energy, where N l is the 
number of fermions in state \i) with Zeeman energy e % z 
determined by the magnetic moments (i B and the mag- 
netic field B. The Zeeman energy term can also be ex- 
pressed as —Hi(N 1 - N 2 ) - H 2 {N 2 - N 3 ) + Ne, where 
the unequally spaced Zeeman splitting in three hyperfine 
levels can be characterized by two independent parame- 
ters Hi — I — €z(fj, B ,B) and H 2 — e%(fi B ,B) — e, with 
I = J2a=i e z(A*B> the average Zeeman energy. 

In general, the scattering lengths depend on spin 
states. However, it is plausible to tune three scatter- 
ing lengths close to each other utilizing the broad Fesh- 
bach resonances 14, 201. Thus the difference in effective 



interaction parameters becomes negligible so that three 
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low spin states may have SU(3) degeneracy. The cou- 
pling constant c/id = —h c/m with interaction strength 
c = —2 /aip determined by the effective ID scattering 
length aiD [2 11 ] - For simplicity, we define a dimensionless 
interaction strength 7 = c/n with density n = N/L. 

The energy eigenspectrum is given in terms of the 
quasimomenta {ki} of the fermions via E = ^ SjLi &| 
which in terms of the function e n {x) = (x + inc/2) /{x — 
inc/2) satisfy the nested BA equations 18, 19j 
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Here j = 1, . . . , N, a = 1, . . . , Mi , to = 1, . . . , M 2 , with 
quantum numbers M x = N 2 + N 3 and M 2 = N 3 . The 
parameters {A Q ,A m } are the rapidities for the internal 
hyperfme spin degrees of freedom. For the irreducible 
representation [3 3 2 2 l Nl ] three-column Young tableau 
encode the numbers of unpaired fermions, bound pairs 
and trions given by Ni = N 1 - N 2 , N 2 = N 2 - N 3 and 
N3 = N 3 , respectively. 

Trions and pairing. In principle, different numbers of 
trions, pairs and unpaired fermions can be selected to 
populate the groundstate by carefully tuning H x and H 2 . 
For a state with arbitrary spin polarization, i.e. N l with 
i = 1,2,3 arbitrary, there are [19( (i) N 3 spin-neutral 
trions in the quasimomentum k space accompanied by 
7V3 spin bound states in the A-parameter space and N3 
real roots in the A-parameter space, (ii) N 2 BCS bound 
pairs in k space accompanied by N 2 real roots in A space 
and (iii) N\ unpaired fermions in k space. With the above 
configuration, we solve the BA equations ((2]) in the strong 
coupling regime L\c\ 3> 1 to give the quasimomenta for 
trions, BCS pairs and unpaired fermions (see Fig. [!}, 
from which the energy is given by 
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Here n a = N a /L with a = 1,2,3 is the density for un- 
paired fermions, pairs and trions, respectively. This state 
can be viewed as a mixture of trionic fermions, hard-core 
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FIG. 1: Schematic configuration of Bethe ansatz quasimo- 
menta k, spin momenta A and A in the complex plane (TV = 29 
with Ni = 6, N2 = 4 and JV3 = 5) at zero temperature. For 
strongly attractive interaction, the unpaired and paired quasi- 
momenta can penetrate into the central region occupied by 
tightly bound trions. 



Eq. ((3]). The root patterns obtained from ([2|) reveal an 
important signature, namely unpaired fermions couple to 
two different kinds - trionic and BCS pairing - of FL. 

Thermodynamic Bethe Ansatz. In the thermodynamic 
limit, i.e. L,N — * 00 with N/L finite, the grand par- 
tition function is Z — tr{e~ u / T ) = e~ G / T , where the 
Gibbs free energy is G = E + Ez — /J,N — TS in terms of 
the Zeeman energy Ez, chemical potential [i and entropy 
S 13, 23, 24 1 . For finite temperatures, besides complex 
BA roots for trions, BCS pairs and real roots for un- 
paired fermions, the quasimomenta {A Q , A m } form com- 
plex strings. The Gibbs free energy can be given in terms 
of the densities of particles and holes for trions, bound 
pairs, unpaired fermions as well as spin degrees of free- 
dom, which are determined from BA l[2|). Thus the true 
physical state is determined by the minimization of the 
Gibbs free energy with respect to these densities, which 
gives rise to a set of coupled nonlinear integral equations 
- the TB A equations [25| ■ 

Quantum phase transitions in the model may be ana- 
lyzed via the dressed energy equations, 
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which follow from the TBA equations in the limit T — > 0. 
Here the function a^ix) = 77- / ■ SM \ -i and are the 

= f + ® a a Jr. - 7/V( a )f 



dressed energies, ctj * e^(x) = f_n" aj(x — y)e^ a \y)dy 



bosons and unpaired fermions, which behave essentially is the convolution. The negative part of the dressed en- 
like particles with different statistical signatures 
The BCS pair binding energy eb = h 2 c 2 /(4m) and the 
binding energy e t = h 2 c 2 /m for a trion can be read off 



ntiaiiy 
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ergies e^(x) for x < \Q a \ correspond to the occupied 
states in the Fermi seas of trions, bound pairs and un- 
paired fermions, with the positive part of correspond- 
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ing to the unoccupied states. The integration boundaries 
Q a characterize the "Fermi surfaces" at e^ a ^(±Q a ) = 0. 
The zero-temperature Gibbs free energy per unit length 
is given by G — Y^,=i S-q" ^ a \x)dx. The chemical 
potential and magnetization per length are determined 
by Hi, H 2 , gm and n through the relations 
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(5) 



In the absence of analytic solutions of Eq. (|U) , we ob- 
tain an exact expansion in the strong coupling regime 
7> 1. Solving the dressed energy equations (g]) by it- 
eration among the relations © and e a (±Q a ) = with 
a = 1,2,3, gives the effective chemical potentials 
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in units of h it /2m. Here we denote /i* = fi + et/S 
for trions, = /U + eb/2 + H 2 /2 for bound pairs and 
fx u = /-i + Hi for unpaired fermions. These results give 
rise to a full characterization of three Fermi surfaces. It 
is important to note [25| that the energy for arbitrary 
population imbalances can be obtained from E/L = /m+ 
G + niHi + n 2 H 2l which coincides with © derived from 
the discrete BA ©. This indicates that the trions and 
BCS bound pairs are possibly the true physical states, 
i.e., the BA roots comprise the equilibrium states in the 
thermodynamic limit. 
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FIG. 2: Phase diagram determined by the energy transfer 
relations (J7J) with chemical potentials © with |c| = 10 and 
n — 1. (a) and (b) show the polarizations m/n and n^/n vs 
the fields Hi and Hi. The figure reveals a novel trion phase 
C, a pairing phase B, an unpaired phase A and four different 
mixtures of these states. 



Full phase diagram. In the strong coupling regime and 
in absence of Zeeman splitting, i.e., Hi = H 2 = 0, the 
dressed energies and e' 1 ' are always positive, i.e., 
Qi = Q2 = 0. Thus trions form a singlet groundstate. 
However, the Zeeman splitting can lift the SU (3) degen- 
eracy and drive the system into different phases. Break- 
ing a trion state requires a spin excitation energy to di- 
minish an energy gap. From Eq. (j4|), the energy transfer 
relations among the binding energy, the Zeeman energy 
and the variation of chemical potentials between different 
Fermi seas are given by 
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These equations determine the full phase diagram and 
the critical fields triggered by the Zeeman splitting Hi 
and H 2 . A similar energy transfer relation was iden- 
tified in exp eriment for ID polarized two-component 
fermions [20(. These relations hold for arbitrary inter- 
action strength. However, for 7 ^> 1, the chemical po- 
tentials are given by Eq. ©. Figure [5] shows the po- 
larization, which clearly indicates novel magnetism, new 
quantum phases, multicritical points and phase transi- 
tions in terms of Zeeman splitting. 

The groundstate energy vs Zeeman splitting parame- 
ters Hi and H 2 can be evaluated from Eq. ([3]) with the 
densities ni and n 2 determined from |7|). Figure [3] shows 
the energy surface for all possible phases shown in Figure 
[2] This figure demonstrates the interplay between differ- 
ent physical groundstates. We see a mixture of trions, 
BCS pairs and unpaired fermions (A + B + C) populates 
the groundstate for certain values of Hi and H 2 . There 
are six different phase transitions across the A + B + C 
boundaries in the Hi~H 2 plane. All phase transitions 
between two phases are second order and reveal a univer- 
sality class of linear-field-dependent magnetization with 
finitely divergent susceptibility in the vicinities of criti- 
cal fields. Our analytical results © do not support the 
square-root-field-dependent behaviour of magnetization 
argued in [23j|, but do agree with results for the attrac- 
tive Hubbard model [2rj |. 

For small Hi (i.e., small splitting between the two 
lower levels), a smooth phase transition from a trionic 
state into a mixture of trions and pairs occurs as H 2 ex- 
ceeds the lower critical value H 2 X (see Fig. [3]). When 
H 2 is greater than the upper critical value H 2 2 , the 
atoms within the two lower states form a pure pair- 
ing phase with SU(2) symmetry and the highest level 
remains unpopulated. In this pure pairing phase, the 
three-level system is reduced to a two-level one. Tri- 
ons and BCS pairs coexist when H 2 X < H 2 < H 2 2 . 

The critical fields H? « ^ (S£ - |£(1 + and 

H f * (¥ + x(! + 3?r)) are uni( l uel y deter " 

mined by the second equation in ijTj). The polarization 
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FIG. 3: Groundstate energy vs Zeeman splitting for |cj = 10 
and n = 1. In the vicinity of the multicritical points and the 
phase boundaries, the energy surface varies continuously with 
Hi and Hz. The insets indicate configurations for the cases 
of small Hi and small H^,- The dashed lines in the insets 
indicate the average Zeeman energy e. 



curve nijn indicates that the phase transitions in the 
vicinities of the critical lines iif 1 and ii| 2 are of second 
order. In addition, the phase transitions B — > A+B — > A 
induced by increasing H\ are reminiscent of those in 
the two-component systems 



26j. We see clearly that 
equally spaced Zeeman splitting H\ — Hi does not favour 
spin-dependent charge states in ID multi-component at- 
tractive fermions. 

For small H2 (i.e., small splitting between the two 
higher levels), the pairing phase is not favored and the 
trions are broken into unpaired fermions in the lowest 
level (see Fig. [3]). Using the first relation in ([7]), we see 
that the trionic state with zero polarization m/n = 
forms the groundstate when the field H < Hf. Here 

ttcI ^ h 2 n 2 ( 2 7 2 7T 2 /-■ , 4 1 

field which makes the excitation gapless. 



is the lowest critical 
If the low- 
est level is widely separated from the two higher levels, 

Hl > H f « ^(22! + 7r 2 (1 __4_ } y all trions 
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are broken and the state becomes a normal FL where 
all atoms occupy the lowest level. Trions and unpaired 
fermions coexist in the intermediate region iff 1 < H 1 < 
Hf. 

For trapped systems, these exotic phases and their seg- 
ments in real space can be predicted with the help of 
the local density approximation. For example, for small 
Hi, a mixture of trions and BCS pairs lies in the trap 
centre and the fully trionic phase (or the fully pairing 
phase) sits in the two outer wings for H 2 < (or 

bh 2 n 2 2 

H2 > — )■ I n contrast, for small H2, a mixture of 
trions and unpaired fermions lies in the centre and the 
fully trionic phase (or the fully unpaired phase) sits in 
the two outer wings for H < e t /3 (or H > e t /3). 

To conclude, our analytic BA and TBA results for the 



critical fields, quantum phase transitions and full phase 
diagrams reveal the nature of spin-neutral trions and 
BCS pairing in three-component ultracold fermions in ex- 
ternal fields. Particularly, we found transitions between 
trionic and BCS pairing phases and between trionic and 
normal FL phases may occur for certain types of Zee- 
man splitting. These exotic phases should stimulate fur- 
ther experimental interest in multi-component ultracold 
Fermi gases with mismatched Fermi surfaces. 
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